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Swift  foxes  ( Vulpes  velox )  are  an  endemic  mesocarnivore  of  North  America  subject  to  resource  and  predation-based  pressures. 
While  swiff  fox  demographics  have  been  documented,  there  is  little  information  on  the  importance  of  top-down  versus  bottom-up 
pressures  or  the  effect  of  landscape  heterogeneity.  Using  a  consumable  resource-based  ideal  free  distribution  model  as  a  conceptual 
framework,  we  isolated  the  effects  of  resource-based  habitat  selection  on  fox  population  ecology.  We  hypothesized  if  swift  fox 
ecology  is  predominantly  resource  dependant,  distribution,  survival,  and  space  use  would  match  predictions  made  under  ideal 
free  distribution  theory.  We  monitored  survival  and  home  range  use  of  47  swiff  foxes  in  southeastern  Colorado  from  2001  to  2004. 
Annual  home  range  size  was  15.4km2,  and  seasonal  home  range  size  was  10.1  km2.  At  the  individual  level,  annual  home  range 
size  was  unrelated  to  survival.  Estimates  of  fox  density  ranged  from  0.03  to  0.18  foxes/km2.  Seasonal  survival  rates  were  0.73  and 
1.0  and  did  not  differ  seasonally.  Foxes  conformed  to  the  predictions  of  the  ideal  free  distribution  model  during  winter,  indicating 
foxes  are  food  stressed  and  their  behavior  governed  by  resource  acquisition.  During  the  rest  of  the  year,  behavior  was  not  resource 
driven  and  was  governed  by  security  from  intraguild  predation. 


1.  Introduction 

Swift  foxes  ( Vulpes  velox )  are  a  mesocarnivore  endemic  to 
the  Great  Plains  region  of  North  America.  Historically,  swiff 
foxes  occupied  the  extensive  shortgrass  prairie  regions  from 
central  Canada  into  New  Mexico  and  Texas  and  from  the 
Rocky  Mountains  east  into  Iowa  [1,2].  Today,  they  are  found 
in  a  variety  of  landscapes,  ranging  from  shrub  steppe  to 
agricultural  to  ranchland  to  native  shortgrass  prairie  [3,  4]. 
While  they  are  capable  of  exploiting  a  wide  range  of  habitats, 
reported  population  parameters  such  as  density  and  survival 
vary  widely  indicating  variation  in  habitat  quality.  There 
is  little  information  on  landscape  structure  or  disturbance 
thresholds  that  control  swift  fox  abundance  or  lead  to 
exclusion  [5]. 

The  ideal  free  distribution  [6]  has  been  used  to  study  the 
influence  of  predation  on  prey  distribution  [7] .  According  to 


the  tenets  of  the  ideal  free  distribution,  high-quality  habitat 
should  be  occupied  to  a  certain  threshold  at  which  point 
competition  or  social  structure  forces  animals  into  poorer- 
quality  areas.  In  this  case,  survival  and  reproductive  success 
would  be  similar  across  landscapes  though  density  and  home 
range  size  would  vary.  However,  most  applications  of  the 
ideal  free  distribution  to  predator-prey  interactions  have 
assumed  predation  risk  is  constant  between  habitats  [8,  9]. 
Predation  risk  is  rarely  constant  across  habitats  and  patches 
[10],  and  extensive  work  has  indicated  that  under  the  threat 
of  predation,  prey  species  preferentially  forage  in  “safer” 
habitats  [10-13].  The  inclusion  of  variable  predation  risk 
into  ideal  free  distribution  models  has  indicated  that  prey 
density  is  more  strongly  influenced  by  the  “riskiness”  of  the 
habitat  than  by  resource  availability  [8]. 

Swiff  fox  population  parameters  in  different  landscapes 
have  been  well  documented  [14];  however,  there  is  little 
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information  available  addressing  the  differences  between 
reported  values.  In  order  to  begin  addressing  the  effects 
of  landscape  structure  on  swiff  fox  distribution  and  demo¬ 
graphics,  we  explored  the  effects  of  landscape  heterogeneity 
on  population  parameters  likely  to  reflect  habitat  quality, 
such  as  population  density,  home  range  size,  and  survival. 
For  this  paper,  we  assessed  the  relative  quality  of  different 
landscapes  as  fox  habitat  based  on  vegetation  structure  and 
prey  base  using  the  ideal  free  distribution  theory  as  a  concep¬ 
tual  framework.  While  the  ideal  free  distribution  is  limited 
due  to  its  resource-based  approach,  it  provides  a  conceptual 
starting  point  for  investigating  the  relative  importance  of 
resource  availability  on  a  native  mesocarnivore. 

Under  the  ideal  free  distribution  theory  we  made  the 
following  series  of  predictions.  (1)  At  the  regional  level, 
estimates  of  swift  fox  density  will  vary  across  landscapes 
while  survival  rates  will  not.  (2)  At  the  individual  level, 
swiff  fox  survival  rates  will  be  unrelated  to  home  range  size. 
(3)  Swift  fox  home  range  size  will  be  negatively  related  to 
rodent  prey  base.  (4)  Swiff  fox  home  range  size  will  be  related 
to  vegetative  structural  elements  that  enhance  rodent  prey 
availability.  The  strength  of  these  relationships  or  variation 
from  these  predictions  can  indicate  the  relative  importance 
of  resource-based  influences  on  swift  fox  ecology.  To  address 
these  predictions,  we  evaluated  swiff  fox  spatial  ecology  at 
4  spatial-temporal  scales  in  southeastern  Colorado  between 
2001  and  2004.  Scales  included  combinations  of  individual 
versus  population  parameters  and  seasonal  versus  annual 
responses.  By  comparing  field  data  to  the  predictions  made 
under  the  ideal  free  distribution  model,  we  were  able  to 
investigate  the  relative  influence  of  top-down  (predation) 
versus  bottom-up  (resource/vegetation)  factors  on  swift  fox 
distribution  and  range  use.  The  study  was  conducted  on  and 
around  the  US  Army  Pinon  Canyon  Maneuver  Site  where 
recent,  dramatic  shifts  in  the  disturbance  regime  due  to 
military  purchase,  as  well  as  discrete  boundaries,  resulted  in  a 
locally  heterogeneous  environment.  The  spatial  scale  of  our 
research  was  small  enough  to  assume  that  prior  to  military 
ownership,  the  landscape  had  been  fairly  homogonous  and 
swiff  fox  population  parameters  were  spatially  consistent. 
The  abrupt  shift  in  ownership,  the  discrete  boundaries,  and 
the  patterns  of  land  use  in  the  area  coalesced  into  a  natural 
experiment  on  the  effects  of  landscape  structure  on  a  native 
mesocarnivore. 


2.  Materials  and  Methods 

2.1.  Study  Site.  We  conducted  research  on  and  around  the 
1,040  km2  Pinon  Canyon  Maneuver  Site  (PCMS)  located 
in  Las  Animas  County,  Colorado  (Figure  1).  North  of  the 
PCMS  the  study  area  extended  into  the  US  Forest  Service, 
Comanche  National  Grassland.  The  study  area  also  extended 
southward  onto  private  ranchlands.  The  region  was  classified 
as  semiarid  grassland  steppe,  with  approximately  60%  of  the 
PCMS  categorized  as  shortgrass  prairie  dominated  by  blue 
grama  ( Bouteloua  gracilis),  western  wheatgrass  ( Agropyron 
smithii),  and  galleta  ( Hilaria  jamesii)  [15].  Shrublands 
interspersed  throughout  the  study  area  included  four- winged 


saltbrush  ( Atriplex  canescens)  and  greasewood  ( Sarcobatus 
vermiculatus),  as  well  as  prickly  pear  cactus  (Opuntia  phaea- 
cantha ),  tree  cholla  ( Cylindropuntia  imbricata ),  and  yucca 
( Yucca  glauca).  The  remaining  landscape  was  dominated 
by  pinyon-juniper  woodland  ( Pinus  edulis  and  Juniperus 
monosperma ).  Elevation  varied  between  1,310  and  1,740  m, 
average  temperatures  ranged  from  1°C  in  January  to  23 °C 
in  July,  and  precipitation  averaged  30  cm  but  can  fluctuate 
widely  [16].  Monthly  precipitation  was  the  highest  in  July 
with  an  average  of  4.3  cm  of  rain,  though  the  35%  of  the 
annual  precipitation  that  fell  during  the  cool-season  (March- 
May)  had  a  proportionally  greater  impact  on  productivity 

[17]. 

The  maneuver  site  was  purchased  by  the  US  Army  in 
1982  to  facilitate  large-scale  mechanized  infantry  training 
activities.  Prior  to  this,  both  the  maneuver  site  and  the 
surrounding  area  had  been  grazed  by  domestic  livestock 
since  the  late  1800s  [18].  Livestock  were  excluded  in  1983 
and  training  began  in  1985  [17].  On  the  maneuver  site, 
training  intensity  varied  widely  simply  due  to  the  logistics 
of  maneuvering  large,  mechanized  infantry  units.  As  a 
result,  remnant  areas  subject  to  very  little  disturbance  were 
scattered  throughout  the  base.  At  the  same  time,  areas 
surrounding  the  base  continued  to  be  grazed  by  domestic 
livestock.  This  natural  heterogeneity  provided  the  basis 
for  our  identification  of  areas  with  variable  vegetation 
structure. 

2.2.  Study  Design.  We  established  6  study  sites  in  areas  sub¬ 
jected  to  3  land  use  regimes:  livestock  grazing,  mechanized 
military  training,  and  unused.  Unused  sites  were  peripheral 
portions  of  the  military  base  that  were  not  used  for  military 
training  yet  were  subject  to  fire  suppression  and  livestock 
exclusion.  Sites  were  named  according  to  local  landmarks 
or  historical  owners:  Private  (PRV),  Biernacki’s  (BTS), 
Pronghorn  (PRN),  Red  Rocks  (RRK),  Bent  (BNT),  and 
Comanche  (COM).  Data  collection  was  based  on  3  seasons: 
winter  (breeding/gestation:  15  December-14  April),  summer 
(pup-rearing:  15  April-14  August),  and  fall  (dispersal:  15 
August-14  December)  [18-20].  Between  2001  and  2004,  we 
conducted  fox  trapping  surveys,  radio  telemetry  of  collared 
foxes,  small  mammal  trapping,  and  vegetation  surveys  on 
each  site  during  each  season  every  year. 

2.3.  Swift  Fox  Capture  and  Radiotelemetry.  We  captured 
swift  foxes  using  double-door  box  traps  (Tomahawk  Live 
Trap  Company,  Tomahawk  WI)  baited  with  chicken  [20, 
21].  Traps  were  placed  500  m  apart  along  a  10  km  transect 
bisecting  each  site  resulting  in  21  traps  per  site.  Each  trap 
was  oriented  and  covered  with  brush  to  provide  protection 
from  exposure.  Traps  were  set  in  the  late  afternoon,  checked 
early  the  following  morning,  and  left  closed  throughout  the 
day.  Each  transect  was  trapped  for  4  consecutive  nights  3 
times  per  year.  For  recollaring  or  targeting  animals,  a  trap- 
enclosure  system  was  used  at  den  sites  [22].  Captured  foxes 
were  handled  without  anesthesia  and  weighed,  sexed,  and 
aged  through  tooth  wear  (adult,  juvenile).  Foxes  were  con¬ 
sidered  juvenile  until  the  pup-rearing  season  following  their 
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( 1)  Private  (grazed)  (4)  Red  rocks  (unused) 

(2)  Biernacki’s  (military  activity)  (5)  Bent  (unused) 

(3)  Pronghorn  (military  activity)  (6)  Comanche  (grazed) 


Figure  1:  Six  study  sites  on  and  around  the  Pinon  Canyon  Maneuver 
as  well  as  the  associated  dominant  land  use  for  that  site. 


birth  (15  April).  Foxes  were  ear- tagged  and  collared  with  30- 
50  g  radio  transmitters  (Advanced  Telemetry  Systems,  Isanti, 
MN). 

We  located  foxes  a  minimum  of  3  times  per  week,  twice 
during  nighttime  hours  when  animals  were  actively  hunting, 
and  once  during  daylight  hours  to  locate  den  sites.  Locations 
were  considered  independent  when  separated  by  at  least  4 
hours  [23],  more  than  sufficient  time  for  a  swift  fox  to  cross 
its  home  range.  Nocturnal  locations  were  estimated  using 
triangulation  of  2-3  bearings  within  5  minutes  and  separated 
by  at  least  40°.  Triangulation  was  done  using  Program 
Locate  II  (Pacer,  Truro,  Nova  Scotia).  Diurnal  locations  were 
collected  visually  by  approaching  the  animal  until  either  a 
den  could  be  identified  or  the  animal  was  seen. 


2.4.  Swift  Fox  Population  Parameters.  Swiff  fox  home  ranges 
were  calculated  at  4  spatiotemporal  scales;  annual  and  sea¬ 
sonal  ranges  at  both  site  and  individual  levels.  All  calculated 
ranges  were  based  on  a  minimum  of  30  locations.  Previous 
area-observation  analyses  [18]  indicated  31  locations  were 
the  minimum  needed  to  estimate  an  individual,  unbiased 
swiff  fox  home  range.  All  ranges  used  in  analyses  were 
based  on  the  95%  fixed  kernel  isopleth  [24]  using  a  least- 
squares  cross-validation  smoothing  parameter.  Least-squares 
cross-validation  allows  the  smoothing  parameter  ( h )  to 
be  chosen  that  minimizes  the  squared  distance  between 
the  fitted  surface  and  the  target  surface  by  varying  the 
bandwidth  [25].  Core  use  areas  (50%  isopleth)  were  also 
calculated  for  comparison  between  sites.  All  ranges  were  cal¬ 
culated  using  the  Animal  Movements  extension  of  Arc  View 
GIS  3.3  [26].  In  this  environment  there  are  no  “habitat 
patches”;  quality  is  a  gradient  blending  from  suitable  into 


Site  in  southeastern  Colorado.  Locations  of  6  transects  are  indicated, 


unsuitable  and  all  home  ranges  contain  elements  of  each. 
Therefore,  kernel  range  estimators  are  more  appropriate 
than  minimum  convex  polygon  methods  due  to  the  fact 
that  there  are  areas  within  ranges  that  are  avoided  by 
residents. 

Site-level  survival  rates,  both  seasonal  and  annual,  were 
estimated  using  the  known  fate  model  in  Program  MARK 
[27].  Due  to  our  emphasis  on  established  breeding  territo¬ 
ries,  only  adult  animals  were  used.  Individuals  not  located 
during  a  season  were  censored  for  that  season.  Individual 
fox  survival  was  calculated  as  the  proportion  of  either 
season  or  year  that  the  fox  survived  [20].  Proportions  were 
weighted  by  the  number  of  days  they  could  be  monitored 
after  capture.  Periods  when  the  animal  was  not  found  were 
censored  due  to  our  inability  to  determine  whether  the 
animal  was  present  and  missed  or  had  temporarily  left  the 
area. 

We  estimated  the  number  of  foxes  in  each  site  using  the 
robust  model  in  Program  MARK.  Seasonal  survival  estimates 
were  taken  from  a  telemetry-based  known  fate  model  due  to 
the  greater  precision.  Dispersal  rates  were  estimated  based  on 
the  number  of  radioed  animals  that  abandoned  an  existing 
home  range.  Juvenile  animals  were  considered  dispersers  if 
they  left  the  vicinity  of  their  natal  territory,  whether  or  not 
they  survived  to  establish  a  range.  Population  estimates  were 
converted  into  density  estimates  by  calculating  the  “effective 
trapping  area”  associated  with  each  transect  [18,  28].  The 
radius  of  the  average  seasonal  95%  kernel  home  range  for 
each  site  was  used  to  buffer  the  transect  in  ArcView  GIS 
[18,  28].  The  resulting  polygon  was  considered  the  “effective 
trapping  area”  for  that  site  [18,  28].  Density  estimates  were 
consolidated  into  seasonal  averages  as  well  as  an  overall 
estimate  for  each  site. 
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2.5.  Landscape  Variables.  For  each  of  the  6  study  sites,  we 
evaluated  vegetation  structure  based  on  50  X  70  m  sampling 
grids  randomly  located  <  1  km  of  the  trapping  transect. 
Four  grids  were  sampled  each  season,  and  new  grids  were 
selected  each  subsequent  season.  Each  grid  consisted  of  seven 
50  m  line  transects  oriented  north-south  and  spaced  10  m 
apart.  On  each  line  transect,  vegetation  type  and  height 
were  evaluated  by  dropping  a  measuring  pin  every  1  m 
and  recording  the  type  and  height  of  the  tallest  vegetation 
encountered  [29].  For  each  grid,  point  measurements  were 
combined  into  estimates  of  percent  basal  cover,  percent  bare 
ground,  percent  litter  (dead  material),  and  mean  shrub  and 
grass  height.  Shrub  density  was  calculated  by  counting  all 
woody  plants  >20  cm  high  within  the  grid.  Grid  estimates 
were  combined  into  seasonal  and  annual  averages  for  each 
study  site.  Standard  deviations  of  grid  estimates  for  each 
study  site  were  used  to  represent  the  homogeneity  of 
vegetation  characteristics  across  each  study  site. 

Following  vegetation  sampling,  35  Sherman  live  traps 
were  placed  throughout  the  grid  with  10  m  spacing  to 
evaluate  swift  fox  prey  base.  Analysis  of  659  swift  fox  scats 
collected  on  the  PCMS  between  1997  and  1998  found  small 
mammals  constituted  the  bulk  of  the  foxes’  diet  throughout 
the  year  [19].  Insects  were  an  ephemeral  food  resource  only 
during  the  summer,  but  rodents  were  utilized  year-round 
and  therefore  we  concentrated  prey  sampling  on  the  rodent 
prey  base.  Traps  were  baited  with  equine  sweet  feed  (corn, 
oats,  and  molasses).  Trapping  grids  were  run  for  4  nights, 
checked  and  closed  each  morning,  and  reset  each  afternoon. 
Captured  rodents  were  marked  with  Sharpie  pens  on  the  tail 
and  abdomen  allowing  for  identification  of  recaptures  over 
the  4-day  trapping  period.  Species  richness  was  calculated 
as  the  number  of  species  captures.  Diversity  was  estimated 
using  the  Shannon-Weaver  index  [30].  Prey  biomass  was 
estimated  using  the  species  average  weight  [31]  multiplied 
by  the  number  of  individual  animals  captured. 

The  Pearson  correlation  coefficients  were  used  to  reduce 
the  number  of  variables  prior  to  analysis.  Final  vegetation 
variables  selected  were  mean  grass  height,  mean  shrub 
height,  percent  basal  cover,  shrub  density,  standard  deviation 
of  basal  cover,  standard  deviation  shrub  height,  and  standard 
deviation  shrub  density.  Standard  deviation  variables  were 
included  to  represent  landscape  heterogeneity.  Final  prey 
base  variables  selected  included  capture  rates  of  the  three 
most  common  species  (Ord’s  kangaroo  rat  ( Dipodomys 
ordii),  northern  grasshopper  mouse  ( Onychomys  leuco- 
gaster),  and  deer  mouse  ( Peromyscus  maniculatus)),  total 
captures  of  all  species,  mean  per  capita  biomass,  species 
richness,  and  species  diversity.  The  standard  deviation  of 
richness  and  diversity  was  included  as  predictor  variable 
to  reflect  community  heterogeneity.  Per  capita  biomass 
was  included  because  as  grassland  converts  to  shrubland, 
a  community  shift  may  sometimes  occur  where  relatively 
rare,  large  bodied  rodents  are  replaced  by  smaller-bodied 
species  with  fewer  antipredator  defenses  [32].  All  vegetation 
and  prey  base  variables  were  log-transformed  to  improve 
normality. 

We  used  AN OVA  to  determine  whether  density  or 
survival  varied  between  sites.  Stepwise  multiple  regression 


was  used  to  identify  which  landscape  variables  best  predicted 
fox  home  range  size  (SAS  v.9.2).  Independent  variables  with 
P  <  0. 1  were  included  in  the  regression  models  and  tolerance 
values  were  set  at  0.1  to  control  for  multicollinearity  [33,  34], 

3.  Results 

A  total  of  47  (24  M :  23  F)  individual  animals  survived  long 
enough  and  met  the  minimum  number  of  locations  for 
calculating  a  seasonal  home  range  at  least  once  (Table  1). 
For  these  47  animals,  we  calculated  69  temporally  indepen¬ 
dent  annual  (34  male  and  35  female)  and  136  temporally 
independent  seasonal  (68  male  and  68  female)  home  ranges. 
Overall  the  average  annual  home  range  size  of  swift  foxes 
was  15.4  ±  9.2  (SD)  km2,  and  the  average  seasonal  home 
range  size  was  10.1  ±  6.7  km2.  The  average  core  use  area 
(50%  kernel  isopleth)  for  swift  foxes  was  1.2  ±  1.1  km2.  The 
average  seasonal  home  range  and  core  use  area  for  male  foxes 
were  12.0  ±  6.9  km2  and  1.5  ±  1.3  km2,  respectively.  Female 
home  ranges  and  core  use  areas  were  smaller,  averaging  8. 1  ± 
5.8  km2  and  0.9  ±  0.8  km2,  respectively.  Estimated  seasonal 
home  ranges  in  the  breeding  season  were  14.2  ±  8.2  km2  for 
male  foxes  and  7.3  ±  4.0  km2  for  female  foxes.  During  the 
pup-rearing  season,  average  home  ranges  were  9.5  ±  5.6  km2 
for  male  foxes  and  6.2  ±  3.6  km2  for  female  foxes.  During 
the  dispersal  season,  ranges  averaged  12.7  ±  86.3  km2  and 
11.7  ±  8.0  km2  for  males  and  females,  respectively.  Home 
range  sizes  and  sample  sizes  varied  among  the  six  sites  (Table 
1). 

3.1.  Prediction  1:  At  the  Regional  Level,  Swift  Fox  Densities 
Will  Vary  Across  Landscapes  While  Survival  Will  Not.  Swiff 
fox  population  density  estimates  (Table  2)  differed  between 
sites  ( F  =  7.48,  P  <  0.001).  Estimates  ranged  from  0.03±0.05 
foxes/km2  on  a  grazed  site  (PRV)  to  0.18  ±  0.10  foxes/km2 
on  an  unused  site  (BNT).  Seasonal  survival  rates  ranged 
between  0.73  on  a  military  site  to  1.0  on  a  grazed  site  and 
rates  did  not  differ  significantly  between  seasons  ( F  =  0.56, 
P  =  0.58).  The  Bent  site  was  excluded  from  this  analysis  due 
to  small  sample  size  and  poor  survival  estimates. 

3.2.  Prediction  2:  At  the  Individual  Level,  Swift  Fox  Survival 
Will  Be  Unrelated  to  Home  Range  Size.  At  the  individual  level, 
annual  home  range  size  was  unrelated  to  survival  (R2  =  0.03, 
F  =  2.68,  P  =  0.11).  Individual-level  results  do  reflect 
some  pseudoreplication  as  more  than  one  annual  range  was 
included  for  15  of  47  animals. 

3.3.  Prediction  3:  Swift  Fox  Home  Range  Size  Will  Be  Negatively 
Related  to  Rodent  Prey  Base.  Between  2001  and  2004,  185 
small  mammal  trapping  grids  were  sampled  for  a  total 
of  740  grid  nights.  Individual  species  capture  rates  were 
uncorrelated  with  each  other,  but  they  were  correlated  with 
rodent  community  indices  such  as  diversity  or  richness. 
The  Pearson  correlation  coefficients  for  these  comparisons 
ranged  from  0.07  to  0.77  (Table  3).  These  correlations  stem 
from  the  fact  that  three  species  (Ord’s  kangaroo  rat,  North¬ 
ern  grasshopper  mouse,  and  deer  mouse)  accounted  for  over 
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Table  1:  Seasonal  home  ranges  and  core  use  areas  for  male  and  female  swift  foxes  in  6  study  sites  located  on  and  around  the  Pinon  Canyon 
Maneuver  Site,  Colorado,  2001-2004.  Parentheses  indicate  sample  sizes  of  foxes  that  were  identical  between  home  range  and  core  use.  Total 
ranges  are  presented  as  ±SE. 


Site 

Sex 

Dispersal 

Home  range  (95%) 
Breeding 

Pup -rearing 

Dispersal 

Core  use  (50%) 
Breeding 

Pup -rearing 

PRV 

M 

10.46  (9) 

11.03  (2) 

6.86  (5) 

1.56 

1.12 

0.86 

F 

7.55  (6) 

7.32  (5) 

4.89  (11) 

0.75 

0.64 

0.72 

BTS 

M 

16.28  (5) 

16.69  (10) 

14.04  (7) 

1.99 

2.00 

2.11 

F 

12.34  (4) 

9.03  (7) 

10.39  (6) 

1.39 

1.16 

0.92 

PRN 

M 

15.94  (5) 

14.71  (3) 

11.09  (5) 

1.58 

1.54 

1.29 

F 

17.33  (6) 

6.65  (4) 

6.33  (6) 

1.88 

0.87 

0.60 

RRK 

M 

13.14  (1) 

8.98  (1) 

6.35  (3) 

0.78 

0.56 

0.75 

F 

18.03  (1) 

6.69  (1) 

6.32  (3) 

1.82 

0.62 

0.68 

BNT 

M 

— 

31.79  (l)1 

— 

— 

4.131 

— 

F 

— 

— 

— 

— 

— 

— 

COM 

M 

7.65  (3) 

6.27  (4) 

4.93  (4) 

1.03 

1.19 

0.60 

F 

5.64  (3) 

2.82  (2) 

2.27  (3) 

0.49 

0.40 

0.24 

Total 

M 

12.63  ±2.0 

12.23  ±  1.6 

8.36  ±  1.0 

1.47  ±0.3 

1.53  ±0.3 

1.04  ±0.2 

F 

11.40  ±  1.8 

7.54  ±  1.5 

6.48  ±  0.9 

1.19  ±0.2 

1.01  ±0.4 

0.66  ±0.1 

throughout  the  study,  only  one  animal  remained  and/or  survived  on 
met  the  requirements  for  calculating  a  seasonal  home  range. 

the  Bent  site  throughout  a 

full  season.  This  animal  was  believed  a 

transient,  though  it 

Table  2:  Estimates  of  population  density  (± 

SE)  and  survival  rates  (± 

SE)  for  adult  swift  foxes  on  6  sites  in  southeastern  Colorado,  2001-2004. 

Site 

Density 

Annual  survival 

Breeding 

Seasonal  survival  rates 

Pup  rearing 

Dispersal 

PRV 

0.18  (0.10) 

0.54 

0.81  (0.09) 

0.83  (0.08) 

0.81  (0.09) 

COM 

0.04  (0.05) 

0.92 

0.92  (0.08) 

1.0  (0.0) 

1.0  (0.0) 

BTS 

0.11  (0.08) 

0.54 

0.84  (0.08) 

0.78  (0.10) 

0.82  (0.08) 

PRN 

0.09  (0.06) 

0.54 

0.73  (0.11) 

0.94  (0.06) 

0.79  (0.09) 

RRK 

0.05  (0.03) 

0.50 

0.80  (0.18) 

0.83  (0.15) 

0.75  (0.22) 

BNT1 

0.03  (0.05) 

— 

— 

— 

— 

throughout  the  study,  only  one  animal  remained  and/or  survived  on  the  Bent  site  throughout  a  full  season.  As  a  result  we  were  unable  to  estimate  survival 
rates  for  that  site. 


75%  of  all  captures  and  therefore  dominated  community 
indices.  Despite  these  correlations,  all  variables  were  used  in 
the  analysis  due  to  biological  importance  and  interpretive 
capacity.  As  stated  above,  regression  tolerance  values  were 
adjusted  to  reduce  the  effects  of  these  correlations. 

Across  all  spatial-temporal  scales,  capture  rate  of  North¬ 
ern  grasshopper  mice  was  the  variable  most  consistently 
associated  with  swift  fox  home  range  size  (Table  4).  Deer 
mouse  capture  rate  was  the  second  most  common  variable, 
emerging  once  as  the  primary  and  twice  as  a  secondary 
explanatory  variable.  While  regression  models  were  generally 
statistically  significant,  they  had  low  explanatory  power  with 
R2  <  0.20.  One  exception  was  the  relationship  between  mean 
population-level  breeding  season  home  ranges  and  Northern 
grasshopper  mice  captures  ( R 1  =  0.467,  P  =  0.02).  Inclusion 
of  community  richness  raised  the  R2  to  0.65  but  raised  the  P 
value  to  0.08.  Core  use  areas,  both  annual  and  seasonal,  were 
unrelated  to  prey  base  variables;  controlling  for  sex  had  no 
effect  on  the  analyses. 


3.4.  Prediction  4:  Swift  Fox  Home  Range  Size  Will  Be  Related 
to  Vegetative  Structural  Elements  That  Enhance  Rodent  Prey 
Availability.  Vegetation  structural  measurements  were  taken 
on  the  same  185  grids  described  above.  A  number  of  the 
measurements  were  correlated,  with  the  Pearson  correlation 
values  ranging  from  0.06  to  0.90  (Table  3).  In  particular, 
percent  basal  cover  and  the  standard  deviation  of  shrub 
height  were  highly  correlated  (P  <  0.05)  with  most  other 
variables.  As  described  above,  all  variables  were  included 
in  the  analysis  due  to  biological  importance  and  regression 
tolerance  values  were  adjusted  accordingly. 

Only  one  vegetative  variable  (shrub  density)  was  consis¬ 
tently  related  to  home  range  size  across  all  spatiotemporal 
scales  (Table  4).  In  general  this  relationship  explained  a  low 
proportion  of  the  variance  in  home  range  size  with  the 
exception  of  the  population  level,  breeding  season  analysis 
(R2  =  0.39,  P  =  0.04).  Core  use  areas,  both  annual  and 
seasonal,  were  not  significantly  related  to  landscape  variables. 
Controlling  for  sex  had  no  effect  on  the  analyses. 
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Table  3:  The  Pearson  correlation  coefficients  for  swiff  fox  prey  base  and  vegetation  structure  variables  on  185  sampling  grids  in  southeastern 
Colorado,  2001-2004.  Capture  rates  indicate  the  total  number  of  individuals  captured  for  selected  species.  Total  captures,  biomass,  and 
community  indices  include  all  species.  Values  exceeding  ±  0.34  are  significant  at  the  P  =  0.05  level. 


Prey  base  variables 

Ord’s 

kangaroo  rat 

N.  grasshopper 
mouse 

Deer  mouse 

Total  captures 

Per  capita 
biomass 

Richness 

Diversity 

SD 

richness 

SD 

diversity 

Ord’s  kangaroo  rat 

— 

0.26 

0.46 

0.57 

0.53 

0.55 

0.62 

0.55 

0.73 

N.  grasshopper 
mouse 

— 

0.18 

0.50 

0.37 

0.50 

0.30 

0.48 

-0.18 

Deer  mouse 

— 

0.46 

0.32 

0.44 

0.50 

0.53 

0.41 

Total  captures 

— 

0.59 

0.99 

0.71 

0.72 

0.48 

Per  capita  biomass 

— 

0.58 

0.41 

0.46 

0.32 

Richness 

— 

0.83 

0.73 

0.56 

Diversity 

— 

0.72 

0.82 

SD  richness 

— 

0.67 

SD  diversity 

— 

Vegetation  structure 
variables 

%  basal  area 

Mean  grass 
height 

Mean  shrub 
height 

Shrubs/ 100  m2 

SD  basal 

area 

SD  shrub 
height 

SD  shrubs/ 100  m2 

%  basal  cover 

— 

0.43 

0.46 

0.30 

0.11 

0.41 

0.23 

Mean  grass  height 

— 

0.76 

0.39 

0.31 

0.20 

0.20 

Mean  shrub  height 

— 

0.65 

0.26 

0.88 

0.52 

Shrubs/ 100  m2 

— 

0.06 

0.39 

0.90 

SD  basal  area 

— 

0.26 

0.15 

SD  shrub  height 

— 

0.35 

SD  shrubs/ 100  m2 

— 

Table  4:  Results  of  stepwise  regression  analysis  of  (a)  vegetation  structure  variables  and  (b)  prey  base  variables  against  home  range  size  for 
swiff  foxes  in  southeastern  Colorado,  2001-2004.  Italics  indicate  a  negative  relationship. 


(a) 


Spatial  scale 

Temporal  scale 

Primary  variable 

Partial 

R2 

Secondary  variable 

Partial  R 1  Model  R2 

Model  P 

Annual 

Shrub  density 

0.044 

None 

0.044 

0.085 

Individual 

Dispersal 

SD  (%  basal  area) 

0.121 

None 

0.121 

0.040 

Pup-rearing 

SD  (mean  shrub  ht) 

0.142 

None 

0.142 

0.028 

Breeding 

SD  (mean  shrub  ht) 

0.284 

None 

0.284 

0.002 

Annual 

Shrub  density 

0.172 

Mean  grass  ht 

0.129  0.301 

0.025 

Population 

Dispersal 

(no  variables  entered  the  model) 

Pup-rearing 

(no  variables  entered  the  model) 

Breeding 

Shrub  density 

0.393 

None 

0.393 

0.039 

(b) 


Spatial 

scale 

Temporal  scale 

Primary 

variable 

Partial 

R2 

Secondary  variable 

Partial  R2 

Model  R2 

Model  P 

Annual 

N.  grasshopper  mouse  captures 

0.086 

Deer  mouse  captures 

0.052 

0.137 

0.056 

Individual 

Dispersal 

N.  grasshopper  mouse  captures 

0.143 

Deer  mouse  captures 

0.138 

0.281 

0.025 

Pup-rearing 

Mean  per  capita  biomass 

0.107 

None 

0.107 

0.068 

Breeding 

deer  mouse  captures 

0.181 

None 

0.181 

0.043 

Annual 

N.  grasshopper  mouse  capture 

0.105 

Total  captures 

0.264 

0.369 

0.001 

Population 

Dispersal 

(no  variables  entered  the  model) 

Pup-rearing 

(no  variables  entered  the  model) 

Breeding 

N.  grasshopper  mouse  captures 

0.467 

Community  richness 

0.179 

0.646 

0.078 
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4.  Discussion 

Currently,  there  are  a  number  of  gaps  in  our  understanding 
of  swift  fox  ecology,  including  the  influence  of  landscape 
variation  on  local  density  and  demographics.  Fundamental 
to  this  understanding  is  the  need  to  quantify  the  influence 
of  resource  limitation  versus  predation  pressure  to  better 
understand  the  processes  influencing  swift  fox  density.  While 
they  may  be  capable  of  exploiting  a  wide  range  of  landscapes, 
swift  foxes  appear  to  be  sensitive  to  thresholds  in  vegetation 
structure,  prey  availability,  or  predator  density  beyond  which 
they  are  excluded.  While  the  ideal  free  distribution  model 
is  limited  in  its  applicability  to  mesocarnivores,  it  is  a 
convenient  point  at  which  to  start  investigating  the  relative 
importance  of  bottom-up  forces. 

Our  results  indicated  swift  foxes  conform  to  2  of  the  4 
predictions  we  made  under  the  ideal  free  distribution  theory 
year-round  but  conform  to  the  remaining  2  only  during 
the  winter  breeding  season.  As  the  ideal  free  distribution 
is  a  resource-based  model,  it  follows  that  swift  foxes  are 
only  resource  (food)  limited  during  the  winter  months. 
During  the  remaining  two-thirds  of  the  year  fox  home  range 
placement  reflected  other  pressures,  most  likely  avoidance  of 
intraguild  predation  by  coyotes  [5,  35].  This  pattern  appears 
to  reflect  what  has  been  termed  an  “adaptive  compromise” 
[11];  during  periods  of  food  stress  swift  foxes  may  trade 
security  for  enhanced  resource  availability.  Similar  behavior 
has  been  well  documented  in  aquatic  systems  and  with 
invertebrates,  birds,  and  mammals;  see  [10]  for  a  review. 

The  idea  that  swift  foxes  make  this  adaptive  compromise 
during  winter  months  is  supported  by  two  additional  lines 
of  evidence.  First,  several  studies  have  reported  that  winter 
months  are  the  most  resource  poor  for  swiff  foxes  [5, 14,  36]. 
In  our  study,  68%  of  small  mammal  captures  occurred 
during  summer  months  (April  15-August  15).  Twenty-three 
percent  of  captures  occurred  in  winter  months  and  9%  in  fall 
months.  The  scarcity  of  mammalian  prey  in  the  fall  is  offset 
by  increases  in  insect  abundance;  insects  made  up  nearly 
80%  of  the  diet  of  swiff  foxes  during  fall  months  [  19] .  While 
mammalian  prey  becomes  important  again  during  winter 
months  [  19] ,  a  significant  portion  of  this  could  be  scavenged 
carrion  [37].  These  results  support  the  idea  that  during 
winter  months  swiff  foxes  may  alter  their  behavior  in  search 
of  additional  food.  Second,  both  our  results  and  a  previous 
study  [18]  suggest  that  the  predation  rate  of  adult  swiff  foxes 
by  coyotes  on  the  PCMS  is  the  highest  in  winter  months. 
While  the  results  were  not  statistically  significant,  the  pattern 
is  consistent  over  5  years.  This  supports  the  idea  that  during 
winter  months  foxes  are  more  exposed  to  predation  risk. 

Swift  foxes  have  the  potential  to  persist  in  a  range 
of  environments,  from  agricultural  to  prairie  remnants  to 
ranchland  [3,  4].  However,  survival,  density,  and  home 
range  estimates  vary  widely  [5].  Interpretation  of  our  results 
suggests  that  swift  fox  behavior  is  a  complex  balance  of  sea¬ 
sonal  food  availability  and  avoidance  of  intraguild  predation. 
Earlier  work  indicated  that  swiff  fox  density  was  significantly 
negatively  related  to  grass  height  [35],  supporting  the  idea 
that  swift  foxes  avoid  areas  of  dense  vegetation  where  they  are 
more  susceptible  to  coyote  predation  due  to  visual  obscurity 


[38].  Other  studies  have  indicated  many  of  the  primary 
prey  items  of  swift  foxes  are  found  in  areas  of  increased 
landscape  heterogeneity  [30,  39].  Therefore,  while  low  shrub 
density  areas  may  be  high-quality,  “safe”  habitat  during  most 
of  the  year,  swift  foxes  in  these  areas  may  need  to  expand 
their  ranges  into  riskier  habitat  during  winter  months  to 
meet  resource  requirements  despite  the  increased  risk  of 
predation.  At  the  same  time,  areas  of  moderate  to  high 
shrub  density  may  provide  for  energetic  requirements  year 
round,  but  the  increased  risk  of  predation  makes  them 
less  desirable  during  spring,  summer,  and  fall  months.  Our 
results  indicate  that  the  highest-quality  swift  fox  habitats  are 
heterogeneous  landscapes  where  resident  animals  can  spend 
the  majority  of  the  year  in  areas  of  low  risk  but  have  access 
to  riskier,  resource-rich  patches  during  winter  months.  Over 
the  past  century,  reductions  in  the  complexity  of  grassland 
disturbance  regimes  due  to  fire  suppression  and  livestock 
husbandry  have  led  to  a  homogenization  of  the  landscape, 
reducing  swift  foxes’  ability  to  effectively  balance  resource 
needs  and  predation  risk. 
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